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Abstract 

Appealing to classical methods of order reduction, we reduce the Lifshitz system to 
a second order differential equation. We demonstrate its equivalence to well known gauge- 
invariant results. For a radiation dominated universe we express the metric and density 
corrections in their exact forms and discuss their acoustic character. 

PACS numbers: 98.80 Hv 
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1 Introduction 



The density perturbations affect the microwave background temperature. The theory of gravi- 
tational instability describes how these inhomogeneities propagate throughout the radiational 
era, and foresee the temperature image they "paint" on the last scattering surface. Classical 
perturbation theory formulated half a century ago by Lifshitz and Khalatnikov ]l|, [|, [3| has 
nowadays been replaced by more appropriate gauge- invariant descriptions |^]- [P1J] . These for- 
malisms introduce some new measures of inhomogeneity. They do not appeal to the metric 
tensor, so they easily avoid spurious perturbations arising from an inappropriate choice of the 
equal time hypersurfaces. They guarantee that the space structures they describe are real 
physical objects. 

On the other hand, the interpretation of the microwave background temperature fluctua- 
tions [12] is based on the Sach- Wolfe effect, where the metric corrections play a key role ]|l3[l . 
Therefore, data obtained from COBE is mostly referred to the classical concepts of Lifshitz 
and Khalatnikov, and only in a minor part to gauge invariant measures, which are more precise 
but difficult to observe Both theories in their original formulations differ essentially. Lif- 
shitz theory provides the two parameter family of increasing solutions for the density contrast 
(|0] formula (115.19)), while all the gauge-invariant approaches foresee in concert only a single 
growing density mode. Thus the interpretation of the microwave temperature map as the initial 
data for cosmic structure formation is fairly ambiguous. 

In this paper we attempt to reconcile both types of theories. We appeal to simple and 
classical methods of order reduction of differential equations HlBl . By use of these techniques we 
remove the pure-gauge perturbations from Lifshitz theory in the radiation dominated universe. 
In consequence we reduce the Lifshitz system to a second order differential equation, exactly the 
same as obtained earlier on the ground of gauge-invariant formalisms. Applying well known 
solutions, we express corrections to the metric tensor, the density contrast and the peculiar 
velocity in exact form. We show that in the early universe, scalar perturbations of any length- 
scale form acoustic waves propagating with the velocity 1 / y/3. 



2 Order reduction 

Relativistic perturbations of a Friedman universe, described in synchronous coordinates [[3], 
0, |3| form a system of two second order differential equations with variable coefficients. In 



contrast, the similar Newtonian problem is expressed by only one second-order equation |16 
|TT| , p~8|] . Obviously, the two additional degrees of freedom appearing in the relativistic case must 
correspond to pure coordinate transformations (gauge freedom) 0, and should be removed from 
the theory. 

Removing pure gauge modes we reduce the Lifshitz equations with pressure p = p/3 to 
Bessel equation. The procedure is as follows: 1) we raise the equations order to fourth, in 
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order to separate the /i n (r]) and X n {v) coefficients, and then 2) we reduce the order of each of 
the separated equations back by eliminating gauge degrees of freedom. The resulting equations 
have exact solutions in the form of Hankel functions H 3 / 2 and their integrals. 

In the synchronous system of reference, the metric corrections h^ u (//, u = 1,2,3) to the 
homogeneous and isotropic, spatially flat universe fulfill the partial differential equations || 
(87rG = c=1) 
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These equations are usually solved by means of the Fourier transform 
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The Fourier transform (^j is defined for absolute integrable functions (the case of least interest 
for cosmology), for nonintegrable functions in the framework of distribution theory, or can be 
understood as a stochastic integral if the initial conditions are given at random [|1^, |2(J . When 
the barotropic fluid (p/p = 5p/5p = w = const) is the matter content of the universe, the 
functions A n (r?) and Pnirf) obey ordinary, second order equations 



- n w(X n (rj) +li n (v)) + 2 
-n 2 w(l + 3w)(\ n (r)) + p n (r])) + (2 + 3w) 



a(rj) 



»n(v) + »n(v) = 0, 
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(5) 



where prime denotes differentiation with respect to the conformal time rj and a is the scale factor 
for the background metric tensor. In order to separate the variable X n (v)) we differentiate (H) 
twice and eliminate terms containing p n {v) or its derivatives by help of eq. (0). We obtain the 
fourth order differential equation 
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In the same way one can treat (H) to find the equation for p n {rf) 
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+ (n 2 w + 2(2 + 3^)^j p'M + (4 + 3w)^p£\r)) + P^(v) = 0. (7) 

In the following part of this paper we restrict ourselves to a universe filled with relativistic 
particles, where both w = Po/po = § and M. = p^a 4 are constants of motion, and the scale factor 



a is a linear function of the conformal time a{rf) = yAi/37]. In the flat universe the expansion 
rate 6{rj) = 3a' \rj) / a{rj) 2 and the energy density Po(r]) relate to each other by po(rj) = 9(t]) 2 /3, 
so the equations for X n {jf) an d Pn(v) take fairly legible form, both prior to 

- + p n (rj)) + -K(v) + K(v) = o, (8) 

3 7] 

2 3 

-n 2 (\ n (r]) + n n {rj)) + -p' n (r]) + p'^7]) = 0, (9) 

V 

and after separation 

(C " ?) x ' niv) + {t + ?) An(r?) + ^ A " )(r/) + A - )(r?) = °' (10) 

^PM + ^p'M + IpKv) + P^iv) = 0. (11) 
67] 3 7] 

We start with equation fllPp. The two well known gauge solutions jl] are (with the accuracy to 
multiplicative constants) 

fi(v) = 1, (12) 
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f 2 (v) = -\/M/3 I —d V = -\og( V ). (13) 



We expect to obtain solutions for (|ToD in the form |l5j : 

Uv) = h(v) (J A( V ) dr^j , (14) 



where A(rj) and -8(77) are some auxiliary functions. Inserting (14-|T5|) into (M) we obtain the 
Bessel equation in its canonical form 

(y ~ ^) + B " [r]) = °- (16) 

Equation ( |T6D is already free of gauge modes, as one can see from simple heuristic considerations. 
Let us assume that there exist a third linearly independent solution of equation (0), which 



corresponds to a pure coordinate transformation. Then, the linear space of gauge modes would 
be 3-dimensional, leaving only a single degree of freedom for the real, physical perturbations. 
Such a theory has no proper Newtonian limit. 

Equation ([Uj) is identical to the Sakai equation ([^TJ formula 5.1), the equation for den- 
sity perturbations in orthogonal gauge (|| formula (4.9), || formulae (16-17)), the equation 



for gauge invariant density gradients (M formula (38)) or Laplacians (M formulae (8-9), JTT 



formula (22)) after transforming these equations to their canonical form (see f22|). It is inter- 
esting to note that equation ([16]) is also identical to the propagation equation for gravitational 
waves ^3f (except for gravitational waves moving with the speed of light). This means that the 
solutions to equation ( |I6"D represent waves travelling with the phase velocity 1/V3 (we show 
this explicitly in the next section). This picture also is consistent 1 with the phonon approach 
as the transformation (f)(i]) = B(r))/r) + B'(rj) to the Field-Shepley variable [ fZ6j , P7[ reduces 



(|i~6|) to the harmonic oscillator (J)"(r]) + \4>(r]) = 



3 Solutions 

The general solution for ([16]) is a combination ofQ 

(17) 




and its complex conjugate, with the frequency u = These solutions are proportional 

to Hankel functions i/3/2, but more frequently are presented as a combination of Bessel and 
Neumann functions b = a\J + a^N ||. Performing integrations (fL4} -|l5|) we determine the 
solution for X n (rj) and find the correction [i n {j]) by solving equation (f|) algebraically 

A(o;77) = Ei(-zW), (18) 

%ijjr\ 

( 1 \ e~ iujr > 

fjt(ur]) = 1 + — + Ei(-iujr)). (19) 

V tujrj J %ior\ 

Obviously, equation ( |TT]) is automatically fulfilled. As a result we obtain the metric corrections 
h^y expanded into planar waves with the frequency constant in conformal time 77 and with 
varying amplitude 

1 The procedure we present here may also be treated as a method to reconstruct metric corrections and 
hydrodynamic quantities in their explicit form, out of the Field and Shepley variables. 
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For similar solutions in the gravitational waves theory see E4 
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+ / A(n)^- 1 + — — + e in -*Ei(-iur)) d 3 n + c.c. (20) 

J 3 \ y iut] J lurj J 

The density perturbation and peculiar velocity can be inferred from formulae (8.2-8.3) of [0 
and expressed as 



— = / A(n)u p (\\ • x, uj7])d 3 n + c.c. (21) 
P 

Sv 



J A(n)u v (n- x,cur])d 3 n + c.c. (22) 

where the Fourier modes form travelling waves 

2 / 1 iujr]\ e*^*-^) , , 

u p (n-x,ur,) = -[1 + — + —L]— , (23) 

u ' {n - x ^ = V3( 1+ ^)^r- (24) 

A generic scalar perturbation in the early universe is a superposition of acoustic waves. Its 
amplitude decreases to reach a constant and positive value at late times. This decrease is 
substantial in the low frequency (early times) limit urj 1. Solutions are formally divergent 
at T] = 0, nevertheless evaluating the cosmic structure backward in time beyond its stochastic 
initiation rji has no well defined physical sense. 

The only perturbations, which are regular at tj — 0, and growing near the initial singularity, 
consist of standing waves u(n ■ x, ujrj) + u(— n ■ x, urf) (compare [^T], [3(J or similar effect 
in the gravitational waves theory [0]). They form a one-parameter family in the 2-parameter 
space of all solutions, so they are non-generic. This property has been confirmed by use of other 



techniques in the gauge-invariant theories [29]. In the stochastic approach nongeneric solutions 



are of marginal interest since they contribute with the zero probability measure. 



4 Summary and conclusions 

It is a matter of dispute whether cosmic structure was created solely by gravity forces JTJ or 
initiated by other, non-gravitational phenomena manifesting themselves as stochastic processes 



jl9| , |20f in some early epochs. For the first hypothesis regular and growing solutions are in- 
dispensable, while in the second one the generic perturbations play a key role. In a radiation 
dominated universe these properties exclude each other. 

Lifshitz theory and the gauge-invariant theories differ less than usually expected. Both types 
of theories, when properly written, lead to the same perturbation equation of the wave-equation 
form. Generic scalar perturbations are superpositions of acoustic waves. Solutions depend on 
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the product ni] (equivalently on urj). Everything which concerns early epochs refers also to long 
waves, and vice versa 3 . The perturbation scale does not divide solutions into different classes. 
Perturbations propagate with the same speed which does not depend on the wave vector. 

This confirms the wave nature of scalar perturbations in the radiation dominated universe (an 
important property already pointed out by Lukash [p5|l , but hardly discussed elsewhere) and 
compels one to use the complete metric corrections (|20|) in the Sachs- Wolfe procedure (not only 
the non-generic growing solutions) at the end of radiational era. 

The reduction technique we apply in this paper can be used for other equations of state. 
For p/p = const 7^ 1/3 solutions can be expressed in terms of hypergeometric functions. In 
other cases solutions may not reduce to any known elementary or special functions, although 
the reduced equation (|16D can be always found. 
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Appendix: Lifshitz "synchronous" gauge 

The original Lifshitz approach |l|, ^|, |3| provides solutions which are different from (p!8|-[T9|), 
and also inconsistent with the gauge-invariant theories. To explain these differences in detail, 
we appeal to the complete solution to (|10|-|l"l|) containing both physical and spurious inhomo- 
geneities. All the gauge freedom within synchronous system is limited to the choice of the 
integral constants in (|T5|). Actually each of these "constants" can be defined as an arbitrary 
function of the wave number n (equivalently to). We write them explicitly as A(~n) and <?(n) 
satisfying 

\(uti) = A (77) (A(n) J A( V ) d V - g(n) log(iw)) (25) 
so they are equal to the Fourier coefficients in the integral 



n,,n 



5 \ ( e*( n ' x ~ wr ?) 



3 This is a peculiar property of the spatially flat radiation-filled universe. 
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e fa - x d 3 n + c.c. (27) 



Each coefficient .A(n), (?(n), can be denned independently. The gauge freedom is carried by 
G(n) what follows directly from (|13|) . Also knowing the gauge invariant methods one can a 
posteriori check that A(n) affects the gauge invariant inhomogeneity measures, while G(n) does 
not. Now, the density contrast and the peculiar velocity, as inferred from formulae (8.2-8.3) of 



— = / [A(n)u p (n ■ x, lot]) + Q{n)u p {n ■ x, cur])} d 3 n + c.c. 

P J 

J [A(n)u v (n ■ x, urj) + Q{jx)u v {xi ■ x, urj)\ d 3 n + c.c. 



(28) 
(29) 



consists of the physical modes u p , u v already found in (^3||2|) and the pure gauge modes equal 
to 



U p (ll • X, LU7]) 
u v (n ■ X, LOT]) 



2 1 



3 rfui 2 
i 1 



(30) 
(31) 



We expand integrals (P8|) and (p9 ) in the early times limit (with the accuracy to i] 2 ), to obtain 

(32) 
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2A(n) + g(n) (\. 1 22,, 



e inx (i 3 n + c.c. 
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e inx d 3 n + c.c. 



(33) 



Both physical and gauge perturbations manifest identical singular behaviour at rj — 0. There- 
fore, one cannot distinguish between them solely on the grounds of their asymptotic forms. On 
the other hand, one is able to regularize perturbations by the gauge choice £?(n) = —A(n). 
Then, the equal time hypersurfaces follow the hypersurfaces of equal density at early epochs. 
This gaugeQ has been actually employed by Lifshitz and Khalatnikov |jj [| , where divergent 
terms l/uj 2 rf are cancelled by the exactly opposite pure-gauge corrections^. In consequence, 
perturbations described there form a mixture of both the physical and the gauge modes. 

4 commonly known as the synchronous gauge 

5 This does not refer to the metric correction where ^-divergence is still present. 
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In the Lifshitz gauge, the mode amplitude [u p (n ■ x, iorj)u p {n ■ x, cur])} 1 ^ 2 grow with time, 
therefore, the two independent solutions for the density contrast increase. The same concerns 
the peculiar velocity. In the low our) limit the density contrast and peculiar velocity form the 
two-parameter linear spaces of growing solutions. As a consequence, a generic inhomogeneity 
increases, which is in conflict with the gauge- invariant theories ||, [5|, 
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